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ABSTRACT 


Due the remarkable growth in the size and complexity of airlift operations, there is an 
increased need for planning tools to assist decision makers with issues ranging from selecting 
the number and types of aircraft for an airlift fleet to making informed decisions with respect 
to investing or divesting in overseas air bases. In Fiscal Year (FY) 94 research was initiated 
in the Operations Research Department of the Naval Postgraduate School in response to a 
request from the United States Air Force Studies and Analyses Agency and resulted in the 
development of a high fidelity strategic airlift optimization model called Throughput II. The 
model is formulated as a multi-period, multi-commodity linear programming model for 
determining the maximum on-time throughput of cargo and passengers that can be 
transported with a given fleet or given network, subject to appropriate physical and policy 
constraints. Troop and equipment movement requirements are specified by the Time Phase 
Force Deployment Data (TPFDD). An optimization model that utilizes the full level of detail 
available in a TPFDD would be of intractable size. Moreover, it is not necessary to build a 
model with such a fine level of detail in order to obtain the important insights required to 
assist decision makers. Therefore Throughput II replaces the potentially large set of airfields 
with a smaller set of centroids and schedules aircraft through these aggregated airfields. 
Currently route selection is performed manually, by an expert, who incorporates a variety of 
factors based on his/her experience. In this thesis we develop techniques for selecting a set 
of candidate routes for any deployment scenario without requiring historical data or extensive 
interaction with an expert. An analyst should be concerned about two potentially detrimental 
effects of these preprocessing procedures. First, infeasibility may be introduced by 
aggregation and second, Throughput II may provide suboptimal solutions since we consider 
a limited number of routes. To address these issues, a postprocessing step can be used to 
screen for constraint violations and to perform sensitivity analysis with respect to alternative 


routing options. 
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EXECUTIVE SUMMARY 


The National Defense Authorization for Fiscal Year (FY) 1991 required the 
Department of Defense to conduct a comprehensive study to determine the mobility 
requirements for the United States Armed Forces. The Mobility Requirements Study 
(MRS) examined all aspects of the mobility question, from domestic transportation to 
intertheater to intratheater requirements. The MRS analyzes sea, air and amphibious lift, 
surface transportation and prepositioning requirements to provide Congress with an 
integrated plan for procuring the necessary lift assets and infrastructure for power 
projection in the 21st century. In January 1991, the Force Structure, Resource and 
Assessment Directorate (J8) of the Joint Staff assumed responsibility for the MRS. The 
Mobility Optimization Model (MOM) [Wing ef al., 1991] was developed by J8 with 
assistance from the Naval Postgraduate School (NPS) within a short period of six weeks. 
MOM was developed to determine airlift, sealift and prepositioning requirements. It 
modeled single onload and offload location with a time window for delivery. Even though 
MOM produces reasonable high level results, it is not suitable for answering the mobility 
questions currently sought by USAF Studies and Analyses Agency (USAF/SAA) because 
it does not model airlift assets, infrastructure and the associated constraints in sufficient 


detail. 


A second linear programming (LP) model for strategic airlift was USAF/SAA’s 


Throughput Model (Yost, 1994). This model was developed to determine an optimal 


xill 


airlift fleet mix based on given movement requirements and to compare different airlift 
fleets in terms of cargo throughput capabilities. Throughput proved to be a useful model 
but had some inherent limitations due to its static nature with respect to time. Important 
system details such as delivery time windows for specific units and unit closure dates are 
not incorporated. In FY 94 research was performed at NPS in response to a request from 
the USAF/SAA, and resulted in the development of a higher fidelity strategic airlift model 
called Throughput II (Lim, 1994) which has since been enhanced (Morton, ef al., 1995). In 
this research, the strategic airlift assets optimization problem is formulated as a multi- 
period, multi-commodity network-based linear programming model, with a large number 
of side constraints. The objective of the model is to minimize late deliveries subject to 
physical and policy constraints, such as aircraft utilization limits and airfield handling 


capacities. 


A typical Time Phased Force Deployment Data (TPFDD) document specifies over 
14000 specific movement requirements over 180 days from 120 Aerial Ports of 
Embarkation (APOE) to 90 Aerial Ports of Debarkation (APOD) through a network 
infrastructure with a large number of potential enroute stops. To attempt to build an 
optimization model that incorporates the full level of detail available ina TPFDD would 
result ina mathematical model of intractable size. Moreover, it is not necessary to build a 
model with such a fine level of detail in order to obtain the important insights required to 
assist decision makers. Therefore, Throughput and Throughput II both use data 


aggregated from the TPFDD for input. Both models, for example, replace the potentially 
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large set of airfields with a smaller set of centroid airfields and schedule aircraft through 


these aggregate airfields. 


We develop a mixed integer location-allocation centroid model that may be used as 
a preprocessing step. Aggregation of airfields may lead to infeasibilities in the resulting 
model with respect to ability to fly certain flight legs due to altered inter-airfield distances. 
This is the principle reason that we propose a centroid aggregation model that has the 
primary objective of minimizing the maximum centroid-satellite distance. A pure minimax 
model, will typically have many multiple optimal solutions. Some optima, however, are 
better than others with respect to criteria ignored by the minimax objective (e.g., the 
satellite-centroid for second furthest airfield). Because of this drawback, a secondary 


minisum objective term is included in the centroid selection model. 


Currently, route selection is performed manually, by an expert, who incorporates a 
variety of factors including critical legs for specific aircraft types, aircrew flight time 
restrictions and commercial track avoidance of certain countries. Such a procedure is time 
consuming and may be adequate for flying to destinations with which USAF has extensive 
experience. However, more difficulties may arise when flying to an unfamiliar destination 
We develop heuristic techniques for selecting the optimal routes for any deployment 
scenario without requiring historical data and extensive interaction with an expert. The 
route selection module is implemented in the Pascal programming language. Feasible 
routes are sorted with respect to their lengths for each origin-destination pair. Inefficient 


and infeasible routes, due to constraints such as crew flying hours or crew stages, are then 
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eliminated. This automated process has the advantage that routes may be generated 
without requiring the extensive experience and knowledge of an expert. Nevertheless, 
there are certain system requirements (e.g., country permissions) that are difficult to 
embed in such a program, and expert knowledge may still be required to rule out the 


feasibility of certain routes 


The elimination and preprocessing step may still generate a large number of routes, 
and it may not be computationally practical to include all the candidate routes in the 
Throughput I] model. For this reason, a method is applied to convert the large problem 
into a smaller problem of manageable size. This reduction is achieved by solving the linear 
program over a subset of its columns; the result of which may be suboptimal. Through 
standard LP sensitivity analysis we can test whether a certain route, not in the original 
subset, is attractive. If it is attractive, it can be put in the basis as a new column and is used 
to improve the route generation model. This procedure can be automated and is known as 
column generation, in our case, the route selection model will be the column generator. 
Including new attractive routes and iterating the Throughput II model can be continued 


until a satisfactory solution is achieved. 
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1. INTRODUCTION 


In spite of the sophisticated inventories and impressive capabilities of modern air 
forces, military aviation is still the youngest of all major military operations. Except for a 
few early attempts, such as using balloons and airships as instruments of war, air power 
has been an extremely effective means of achieving military objectives, mainly seins of 
its great versatility, for the last 80 years. Among their many missions, airforces provide 
two key tasks: delivering firepower directly by air-to-surface and air-to-air weapons, and 
rapidly deploying combat units over substantial distances. While the use of military air 
power ts relatively new, air mobility operations are of even more recent vintage. The 
substantial exploitation of military airlift dates only from 1938. Prior to that time, except 
for some primitive transportation efforts on a limited scale via air transport systems, the 
deployment and logistic support of military forces was undertaken mostly via surface 


transport systems. 


There are three fundemental disadvantages of surface deployment, which cannot be 
eliminated by speed and efficiency provided by the newly developing ideas and technology 
of this century. Firstly, surface deployment over substantial distances is slower than air 
transport systems. Secondly, surface deployment is heavily restricted by geographical 
constraints The third, and most important reason, concerns its susceptibility to enemy 


attacks. On the other hand, the airlift has less comparative capacity than the sealift does, 





but because of its speed and range, all the modern armed forces have incorporated airlift 
capability into their transportation systems. 

During World War II, in addition to the effectiveness of fighting, attacking and 
bombing units of air forces, air transport operations proved to be an efficient means of 
projecting forces to strategic locations. This added another aspect to the obvious potential 
of air forces as a direct force, namely the use of air forces as an indirect force. The 
importance of this latter feature of air forces is easily recognized by the success of airlift 
operations which were totally impossible by other means in the same range and time 


frames 


Despite prepositioned fuel, ammunition, and equipment, the magnitude of the 

United Nations/ United States airlift effort during Desert Shield and Desert Storm was 
unprecedented. By 10 March 1991, strategic airlift had moved over 500,000 people and 
540,000 tons of cargo to the theater. At the height of the Desert Shield airlift, Military 
Airlift Command’s (MAC’s) cargo movement averaged 17 million ton-miles per day’. By 
comparison, during the 1973 Arab-Israeli War, US airlift moved 4.4 million ton-miles per 
day. Other historical comparisons include the World War I] “Hump” at 0.9 million ton- 
miles per day, the Berlin Airlift at 1.7 million ton-miles per day, and Operation Just Cause 


at 2.0 million ton-miles per day (Gulf War Air Power Survey, 1993). 





" US strategic airlift requirements and capabilities may be measured in terms of million-ton-miles per 
day, i.c., in multiples of the capacity to move one ton of cargo by air a distance of one mile in one day. 
This standard of measurement, though a useful planning tool and gauge of present and future strategic 
airlifi capabilitics, does not, of course, take into account such real world constraints as aircraft and crew 
exhaustion, weather, availability of airfields, overflight rights, and possible enemy action. 


Due to the remarkable growth in the size and complexity of airlift operations, there 
is increased need for planning tools to assist decision makers with tssues ranging from 
selecting the number and types of aircraft for an airlift fleet to making informed decisions 
with respect to investing or divesting in overseas air bases. Although it must be interpreted 
carefully, this “search for optimality” is a very important theme in Operations Research. 
The field of Operations Research provides tools for attempting to find the best, or optimal, 


solution to mathematical models of the problem under consideration. 


A. BACKGROUND 


The National Defense Authorization for Fiscal Year (FY) 1991 required the 
Department of Defense to conduct a comprehensive study to determine the mobility 
requirements for the United States Armed Forces. The Mobility Requirements Study 
(MRS) examined all aspects of the mobility question, from domestic transportation to 
intertheater to intratheater requirements. The MRS analyzes sea, air and amphibious lift, 
surface transportation and prepositioning requirements to provide Congress with an 
integrated plan for procuring the necessary lift assets and infrastructure for force 
projection in the 21st century, In January 1991, the Force Structure, Resource and 
Assessment Directorate (J8) of the Joint Staff assumed responsibility for the MRS. The 
Mobility Optimization Model (MOM) (Wing e/ al/., 1991) was developed by J8 with 
assistance from the Naval Postgraduate School (NPS) within a short period of six weeks. 
MOM was developed to determine airlift, sealift and prepositioning requirements. For 


airlift, it modeled single onload and offload locations with a time window for delivery. 


Even though MOM produces reasonable high level results it is not suitable for answering 
the mobility questions currently sought by USAF Studies and Analyses Agency 
(USAF/SAA) because it does not model airlift assets, infrastructure and the associated 
constraints in sufficient detail 

A second linear programming (LP) model for strategic airlift was USAF/SAS’s 
Throughput Model (Yost, 1994). This model was developed to determine an optimal 
airlift fleet mix based on given movement requirements and to compare different airlift 
fleets in terms of cargo throughput capabilities. Throughput proved to be a useful model 
but had some inherent limitations due to its static nature with respect to time. Important 
system details such as delivery time windows for specific units and unit closure dates are 
not incorporated. In FY 94, research was performed at NPS in response to a request from 
the USAF/SAA, and culminated with development of a higher fidelity strategic airlift 
model called Throughput II (Lim, 1994). In this research, the strategic airlift assets 
optimization problem is formulated as a multi-period, multi-commodity network-based 
linear programming model, with a large number of side constraints. The objective of the 
model is to minimize late deliveries subject to physical and policy constraints, such as 


aircraft utilization limits and airfield handling capacities. 


B. PROBLEM STATEMENT 
Executing a deployment - delivering units and equipment to theater efficiently and 
effectively - with limited aircraft and infrastructure assets is a formidable task. The 


National Command Authority (NCA), the Joint Chief of Staff, and major field 


commanders have access to the Worldwide Military Command and Control System 
(WWMCCS), which includes a command and control system called Joint Operational 
Planning and Execution System (JOPES). The Operationai Plan (OPLAN) is entered as a 
requirement and JOPES generates a Time Phased Force Deployment List (TPFDL). This 
list called the Time Phased Force Deployment Data (TPFDD) document, and contains 
information pertaining to the force and movement requirements of the OPLAN. With this 
information, units able to supply force requirement can be identified and transportation 
arrangements can be made. A typical TPFDD specifies over 14000 specific movement 
requirements over 180 days from 120 Aerial Ports of Embarkation (APOE) to 90 Aerial 
Ports of Debarkation (APOD) through a network infrastructure with a large number of 
potential enroute stops. To attempt to build an optimization model that utilizes the full 
level of detail available ina TPFDD would result in a mathematical model of intractable 
size. Moreover, it is not necessary to build a model with such a fine level of detail in order 
to obtain the important insights required to assist decision makers. Therefore, Throughput 
and Throughput II both use data aggregated from the TPFDD for input. Both models, for 
example, replace the potentially large set of airfields with a smaller set of centroid airfields 


and schedule aircraft through these aggregate airfields. 


A mixed integer programming model for centroid selection has been developed by 
AFSAA that may be used as a preprocessing step for either model (Yost, 1994). The basic 
idea in this model is to first define candidate centroid airfields based on their passenger and 


cargo flows and then to solve the model that minimizes a weighted sum of distances 


between the centroids and their satellite airfields. A primary shortcoming of such an 
ageregation procedure is that it may lead to infeasible solutions in the subsequent airlift 
optimization model Airfields aggregated into centroids may inappropriately “borrow” 
capacity from one another. As a result, it may be discovered on disaggregation that the 
solution to Throughput II has assigned too many sorties to an airfield. Such an outcome 
may occur, for example, when a relatively small capacity airfield with a large flow 
requirement is aggregated with another airfield that has a large capacity value. A second 
potential drawback in the current method of aggregating airfields is that the centroids are 
selected to minimize the sum of weighted distances to all corresponding satellite airfields, 
and the distances between airfields in Throughput II are determined from the centroid not 
the satellite. As a result, there are cases in which a particular aircraft type may be capable 
of flying a route which originates from the centroid but not from the satellite airfield. 
Conversely, certain routes may be eliminated unnecessarily due to airfield aggregation. As 
a result of these aggregation effects, the fidelity of the associated Throughput IJ model 


may be reduced 


Currently, route selection is performed manually, by an expert, who incorporates a 
variety of factors including minimum payload ferry ranges” for specific aircraft types, 
aircrew flight time restrictions and commercial track avoidance of certain countries. Such 


a procedure is time consuming but may be adequate for flying to destinations with which 





* This is the longest distance which an aircraft can fly with minimum reasonable payload (20 % of its 
maximum capacity) without refueling. 


USAF has extensive experience; however, difficulties may arise when flying to an 


unfamiliar destination. 


In summary, the size of the infrastructure network is currently reduced by a 
preprocessing step that selects airfields via the centroid model and eliminates, by hand, a 
large number of infeasible or inefficient routes. The results of Throughput IT have 
indicated that for some important deployment scenarios the airlift system may be 
infrastructure constrained, in part, due to the aircraft handling capability of certain enroute 
airfields. This makes the process of preselecting centroid airfields as well as delivery and 
recovery routes very important. It is desirable to minimize the effect of capacity 
borrowing. Moreover, it may be possible to improve the system’s throughput by proper 
aircraft-route scheduling. An automated system that can generate routes for given origin- 
destination pairs could dramatically aid the analyst in the selection of deployment and 


recovery routes 


C, METHODOLOGY 


This thesis will consist of two main parts as illustrated in Figure 1.1. The first part 
is a preprocessing procedure that selects a set of centroids and candidate routes that will 
be used as input to Throughput II or any other appropriate optimization model. The 
second part is a postprocessing procedure that seeks to determine whether the 
Throughput I solution may be disaggregated into a valid solution for the original model 


and/or improved by providing Throughput IT with additional routing options. 
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Figure 1.1. Route Selection Model 


1. Preprocessing Modules 


The preprocessing step.consists of two distinct modules; first a model is solved to 
determine the centroids and then a heuristic® generates routes from the origin centroid 
airfields through enroute airfields, to the destination centroid airfields for all relevant 


aircraft types. 


a. Centroid Model 
We develop a location-allocation model with a primary objective to 
minimize the maximum centroid-satellite distance and a secondary objective to minimize 


the sum of the corresponding distances. Aggregation of airfields may lead to infeasibilities 





3“ heuristic is a rule of thumb, strategy, or trick used to improve the efficiency of a system which tries to 


discover the solutions of complex problems.”’(Slagle, 1971). 


in the resulting model, with respect to an aircraft’s ability to fly certain flight legs, due to 
altered inter-airfield distances This is the primary reason that we propose a centroid 
aggregation model that has a primary objective of minimizing the maximum centroid- 
satellite distance. A pure minimax model, typically has many multiple optimal solutions 
Some optima, however, are better than others with respect to criteria ignored by the 
minimax objective (e.g., the satellite-centroid for second furthest airfield). Because of this 
drawback, the secondary minisum objective term is included in the centroid selection 
model 

b. Route Selection Model 

The route selection module is implemented in the Pascal programming 
language. Feasible routes are sorted with respect to their lengths for each origin- 
destination pair. Inefficient and infeasible routes, due to constraints such as crew flying 
hours or crew stages, are then eliminated. This automated process has the advantage that 
routes may be generated without requiring the extensive experience and knowledge of an 
expert. Nevertheless, there are certain system requirements (e.g., country permissions) 
that are difficult to embed in such a program, and expert knowledge may still be required 


to rule out the feasibility of certain routes 


z Postprocessing Module 


The final output of the preprocessing step is a set of candidate routes that can be 
used as input data for Throughput II or any other appropriate model such as the 


Stochastic Airlift Optimization Model (Goggins, 1995). However the preprocessing step 


may still generate a large number of routes, and it may not be computationally practical to 
include all the candidate routes in the Throughput II model. For this reason, a method is 
applied to convert this large problem into a smaller problem of manageable size. This 
reduction is achieved by solving the linear program over a subset of its columns, the result 
may be suboptimal. Through standard LP sensitivity analysis we can test whether a certain 
route, not in the original subset, is attractive. If it is attractive, it can be put in the basis as 
a new column and is used to improve the route generation model. This procedure can be 
automated and is known as column generation; in our case, the route selection model will 
be the column generator. This process of including new attractive routes in the 
Throughput I] model, solving the corresponding LP, and testing for optimality can be 


continued until a solution to the original problem is achieved. 


D. SUMMARY 


In this thesis we develop techniques for selecting the optimal routes for any 
deployment scenario without requiring historical data and extensive interaction with an 
expert. The results of Throughput II have indicated that for some important deployment 
scenarios the airlift system may be infrastructure constrained. This makes the process of 
preselecting centroid airfields very important. In the first part of this thesis we develop an 
efficient procedure for solving the combination minimax-minisum facility location model 


to determine centroid airfields that will help curb detrimental aggregation effects. 
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An analyst should be concerned about two potentially detrimental effects of the 
preprocessing procedure. Either infeasibility may be introduced by aggregation or 
Throughput II may yield a suboptimal solution since only a limited number of routes are 
considered To address these issues, a postprocessing step includes 
i. disaggregation and feasibility tests, 

il, sensitivity analysis with respect to alternate routing options. 
Together, these measures assist Throughput IJ in providing tractible, as well as 


valid results 
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Il. REVIEW OF THROUGHPUT IT MODEL 


The Throughput I] strategic airlift optimization model was developed by Captain 
Lim Teo Weng (Lim, 1994) and has since been enhanced in ongoing research at the Naval 
Postgraduate School (Morton e/ a/., 1995). The mathematical formulation of Throughput 
I] is contained in Appendix A. The Throughput I] model is a multi-period model for 
determining the maximum on-time throughput of cargo and passengers that can be 
transported with a given fleet over a given network, subject to appropriate physical and 
policy constraints. The objective function minimizes the total weighted penalties for late 


deliveries and nondeliveries 


The preprocessing and postprocessing techniques developed in this thesis are 
utilized by the Throughput II optimization model as indicated in Chapter I. For this 


reason, we provide a review of Throughput I] in this chapter 


A. MODEL FEATURES 


The Throughput II model has been designed to handle many of the airlift system’s 
particular features and modes of operation. The model is a strategic airlift model, meaning 
that it considers inter-theater but not intra-theater deliveries. The major features of the 
airlift system currently captured by the model include (Morton e/ a/., 1995): 


e Multiple origins and destinations: The model routes aircraft through multiple origin, 
enroute and destination airfields. 


Flexible routing structure: The air route structure supported by the model includes 
delivery and recovery routes with a variable number of enroute stops (usually between 
zero and three). This gives the model the option of short-range flights with heavier 
loads or long-range flights with lighter loads. For further routing flexibility, the model 
also allows the same aircraft to fly different delivery and recovery routes on the same 
mission 


Aircraft-to-route restrictions: The user may impose aircraft-to-route restrictions; e.g., 
military aircraft may only use military airfields for enroute stops. This particular 
provision arises because the USAF Air Mobility Command (AMC) may call upon 
civilian commercial airliners to augment USAF aircraft in a deployment, under the 
Civil Reserve Airlift Fleet (CRAF) program. The model distinguishes between USAF 
and CRAF aircraft 


Aircraft assets can be added over time. This adds realism to the model because CRAF 
and other aircraft may take time to mobilize and are typically unavailable at the start of 
a deployment 


Delivery time windows: In a deployment, a unit is ready to move on its available-to- 
load date (ALD) and has to arrive at the theater by its required-delivery-date (RDD). 
This aspect of the problem has been incorporated in the model through user-specified 
time windows for each unit. The model treats the time window as “elastic” in that 
cargo may be delivered late, subject to a penalty 


ASSUMPTIONS 


The assumptions used in the model are as follows: 


Airfield capacity is represented by a single aggregate figure, called Maximum-on - 
Ground (MOG). The literal translation of MOG as the maximum number of planes 
that can be simultaneously on the ground at an airfield is misleading, because MOG is 
used to convey more than just the number of parking spaces. In actuality, airfield 
capacity depends on many additional factors such as availability of material handling 
equipment and various ground services capacities. Unfortunately, data are not 
currently available to support a multidimensional MOG modeling enhancement 


Inventoried aircraft at origin and destination airfields are considered not to affect the 
aircraft handling capacity of the airfield. This assumption is not strictly valid since an 
inventoried aircraft takes up parking space even if it is not consuming services. 


e Deterministic ground time: Aircraft turnaround times for onloading and offloading 
cargo and enroute refueling are assumed to be known constants, although they are 
naturally stochastic. This ignores the fact that deviations from the given service time 
can cause congestion on the ground. To offset the optimism of this assumption, an 
efficiency factor is used in the formulation of airfield capacity constraints. 


Cc. CONCEPTUAL MODEL FORMULATION 


The primary decision variables are the number of sorties initiated, and the amount 
of cargo and passengers carried, for each unit, by each aircraft type, via each available 
route, in each time period. Additional variables are defined for the recovery flights, for 
aircraft inventoried at airfields, and for the possibility (at high penalty cost) of not 


delivering required cargo or passengers 


le Objective Function 


The objective function minimizes the total weighted penalties for late deliveries 
and nondeliveries subject to appropriate physical and policy constraints. The penalties are 
weighted according to two factors: the priority of the unit and the degree of lateness. The 
penalty increases with the amount of time late, and non-delivery has the most austere 


penalty 


The anticipated use of the model is for situations when the given airlift resources 
are insufficient for making all the required deliveries on time. On the other hand, if there 
are enough resources for complete on-time delivery, then the model’s secondary objective 


function is to choose a feasible solution that maximizes unused aircraft. The motivation for 


the secondary objective is that if the available aircraft are used as frugally as possible, 
while still meeting the known demands and observing the known constraints, then the 
mobility system will be as well prepared as possible for unplanned breakdowns and 


unforeseen requirements, such as an additional nearly simultaneous regional contingency 


2. Constraints 


The model’s constraints can be grouped into five categories: demand satisfaction, 


aircraft balance, aircraft capacity, aircraft utilization, and airfield handling capacity 


Demand Satisfaction Constraints: The cargo demand constraints attempt to ensure for 
each unit that the correct amount of cargo moves to the required destination within the 
specified time window. The passenger demand constraints do the same for each unit’s 
personnel. The demand constraints have elastic variables for late delivery and non- 
delivery. The optimization will seek to avoid these options if it is possible with the 
available assets, or to minimize them if not. 


Aircraft Balance Constraints: These constraints keep physical count of aircraft by type 
(e.g., C17, C5, C141, etc.) in each time period. They ensure that the aircraft assets are 
used only when they are available. 


Aircraft Capacity Constraints: There are three different kinds of constraints on the 
physical limitations of aircraft -- troop carriage capacity, maximum payload, and cabin 
floor space -- which must be observed at all times. 


Aircraft Utilization Constraints: These constraints ensure that the average flying hours 
consumed per aircraft per day are within AMC’s established utilization rates for each 
aircraft type 


Aircraft Handling Capacity at Airfields: These constraints ensure that the number of 
aircraft routed through each airfield each day is within the airfield’s handling capacity. 


D. LIMITATIONS 


Two limitations of the optimization model are the inability to handle local traffic 


congestion and low fidelity with respect to certain aspects of the airlift system 


1. Time Resolution 
With one-day time periods, the model can route aircraft in a manner that causes 
local congestion. Another limitation is that the model rounds the time (to the nearest day) 


at which an aircraft arrives at enroute and destination airfields. 


2. Aircraft Reliability 
Aircraft in need of repair can have an immediate impact on throughput capability, 


especially when airfield capacity resources are limited. 


3. Deterministic Ground Time 


The aircraft ground times used in the calculation of MOG consumption represent 


the expected times for onload/refuel/offload, resulting in optimistic throughput capability. 


4. Airfield Aggregation 


Throughput II uses data aggregated from the TPFDD for input. It replaces the 
large set of airfields with a smaller set of centroid airfields and schedule aircraft through 


these aggregate airfields 
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5. Limited Route Structure 


Route selection is performed manually, by an expert, who incorprates a variety of 
factors including minimum payload ferry ranges for specific aircraft types, aircrew flight 
time restrictions and commercial track avoidance of certain countries. Such a procedure is 
time consuming but may be adequate for flying to destinations with which USAF has 
extensive experience, however difficulties may arise when flying to an unfamiliar 


destination 


Airfield aggregation and limited route structure are examined in this thesis. 
Limitations involving aircraft reliability and stochastic ground times are examined by 


Goggins (1995) 


Hil. PREPROCESSING MODULES 


The preprocessing step consists of two distinct modules where airfield aggregation 
and route selection are performed. The first module aggregates origin and destination 
airfields into selected centroid airfields. The second module generates the route set which 


will be used in the optimization model. 


A. CENTROID MODEL 


As outlined in the introduction, aggregation of airfields may lead to infeasibilities in 
the solution produced by the resulting strategic airlift model with respect to borrowing 
airfield capacity and with respect to an aircraft’s ability to fly certain flight legs due to 
altered inter-airfield distances. It is primarily for the latter reason that we propose a 
centroid aggregation model that has a primary objective of minimizing the maximum 
centroid-satellite distance. A pure minimax model will typically have many multiple 
optimal solutions, but some optima are usually better than others with respect to criteria 
ignored by the minimax objective (e.g., the satellite-centroid for second furthest airfield). 
Asa result we incorporate a secondary minisum objective term in the centroid selection 
model to resolve multiple optimal solutions. Aggregation is not performed for enroute 
airfields, but is performed for origins and destinations separately. 

The preprocessing step is based on a set of raw data provided by various sources 
such as the TPFDD and airfield infrastructure data. Due to some policy constraints or data 


inaccuricies, there are occasionally some triangle inequality violations for inter-airfield 
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distances. As a data integrity check we ensure that given any two points X = (x, mal 


and Y =(y, >) the great circle distance between X and Y, which we denote by |X - Y' has 


the following properties for all X and Y: 





1. Nonnegativity Ix -¥/>0 
2 Symmetry |X -¥|=|Y- | 
3. Triangle Inequality |X -Y|<|X -Z|+|Z-Y| for any Z 


For the given data set, if the triangle inequality violation occurs, 
|X -¥|>|X-Z|+|Z-¥Y| for any X, ¥, Z 
an algorithm is used to replace this inequality with the following equality 


|x -Y=|x -z|+|z-Y| 


1. Mathematical Formulation 


a. Index Sets 
ié/_: Original airfields, e.g., Travis AFB 
j €J * Candidate airfields that may be centroids 
b. Data 
d, Great circle distances between airfield / and centroid j 
k . Number of origin (or destination) airfields that will be used in 


the Throughput IT model 
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G Decision Variables 


1 of airfield ts assigned to centroid j 
xX, = 
’ 10 otherwise 


| if candidate airfield j is selected asa centroid 
119 otherwise 


a Minimax Model 


Kt. SY, V(i,j)elxJ 
2,0, >h 
jes 
x, =1 Viel (2.1) 


x, €{0,]} Vi, j)elxJ 
y, €{0,1} VieJ 
The integer program in (2.1) is a variation of Schrage (1991) that may be 


converted to a linear integer programming problem by introducing a new decision variable 


zvia 


min 2z 

Sh 22d, +x, V(i,j)elxJ 
x, Sy, Vij) El x J 
Dy, =k @.2) 
yes 


x =] Viel 
pas 


x, € {0,1 V(i,j) EL x J 
y, €{0,1} VieJ 


2) 


This model minimizes the maximum distance between centroids and their 
assigned airfields while ensuring that all airfields are assigned to only one centroid and the 
total number of centroids equals a prespecified value. For reasons detailed above, we use a 
secondary minisum objective to select among multiple optimal solutions. This can be 


stated mathematically, 


z’ =min z+ey> > d,.*%, 


wl yes 


Rh fb. +% Vi, j)elxJ 


yoy 


x, Sy, V(i,j)elxJ 

Sy, =k (2.3) 
ied: 

>> x, =1 Viel 

ye 


x, € {0,1} Vii,j)elxJ 
y, €{0,1} vied 


where € > 0 1s a sufficiently small prespecified constant. Model (2.3) contains a large 
number of binary variables that can quickly result in unreasonably long solution times, 


even for modest values of |/| and |./|. As a result, we utilize a solution technique that 


solves a sequence of a set covering location models to find the optimal solution to (2.3). 
é. Set Covering Model 


The covering problem involves selecting the minimum number of centroids 


so that each of the original airfields is within a specified distance of some centroid (i.e., 
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covered by that centroid) (Francis e¢. al., 1992). The set covering model can be stated by 


introducing the following addition notation, 


(1) Data 
r Maximum allowed distance between a centroid and its assigned 
airfields 
(2) Index Set 


oO = { jid,sr \ set of centroids that may used to cover airfield | 


z= min) y, (2.4) 


yes 


SE. ere Viel 


102, 


y,e{0]} Vie 


The key connection between the set covering model (2.4) and minimax / € -minisum model 
(2.3) is the following. If Leni , the optimal objective value to (2.4) for a specific value of r, 
exceeds k ( the centroid budget in (2.3)) then we know z/,,) > r. Similarly if 27, ,. <k 
then we know z/,,, <7. Our goal, then, is to find the smallest value of r for which (2.4) 


yields an objective of at most &. Such an approach proves to be computationally efficient 


because: 


(i) the number of binary decision variables in (2.4) is |J| while the number of binary 


variables in (2.3) is |/| + |./| + |J| and 
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(ii) empirical evidence has shown that a linear relaxation of (2.4) often provides a “near- 
integer” solution. 
The algorithm for solving (2.3) via a sequence of set covering models as defined in 


(2.4) is shown in Figure 2.1. The main part of the algorithm solves (2.3) with ¢ = 0, anda 


postprocessing step solves an appropriately defined minisum model to resolve multiple 


optimal minimax solutions according to the secondary minisum objective. 
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Figure 2.1 Set Covering Location Algorithm 
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The post-processing minisum model may be stated : 


Index Set 


Cy = {@/)id, <r \ 


$t. %, i, V(i, /) €Q 
y= 
gua 
> x, =! Viel (2.5) 
yd 
y, €{0,1} VieJl 
O<x, <1 V(i,s) EQ 


Note that due to the linearity of the objective the binary constraints on x; can be 
telaxed in (2.5) (Schrage, 1991). This mixed integer program selects a solution from the 
set of optimal solutions to (2.1), because centroid assignments are approximately 
restricted by the set 2. General Algebraic Modeling System (GAMS) (Brooke ez. al, 
1992) code for solving model (2.3) via the algorithm of figure 2.1 is in Appendix B. 


Computational results for this and alternative approaches are discussed in Chapter V. 


B. ROUTE SELECTION MODEL 

The primary goal of strategic airlift planning is to satisfy Commander-in-Chief 
(entral Command (CINCCENT)’s requirements by employing airlift resources effectively. 
To meet the goal, planners have to consider the entire airlift system and its interrelated 


parts. Each enroute stop and destination airfield must have adequate runways, taxiways, 
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ramps and support facilities. Where crew changes are required, a sufficient number of 
qualified and properly rested aircrews must be available. Normal crew duty time 
limitations (unless otherwise specified in the governing OPORD/OPLAN) are 16 hours for 
a basic crew and 24 hours for an augmented crew. 

The route selection model may be viewed as a search process or more specifically 
as a traversal of a tree structure in which each node represents an airfield and each arc 
represents a flight leg between the airfields represented by the nodes. The search process 
must find paths through the tree that connect the origin with the destination and obey 
given physical constraints for each aircraft type such as the appropriate airfield types it 
may utilize, the aircraft minimum payload ferry range, crew stages etc. Because of the 
extremely large number of possible routes that could be taken, we need to bound the 
number of enroute stops. Governing our search is a set of rules that defines the legal 
routes we can follow. These rules include the crew stages, enroute airfields, runway 
limitations in airfields and the type of airfield (e.g., military, civilian, joint), Suppose that 
we want to find some routes from node A to node G through a network of airfields, such 


as shown in Figure 2.2 

















Figure 2.2 A Basic Route Selection Problem. 
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Finding a route involves six kinds of efforts: 

e Each leg of the flight should be less than the minimum payload ferry range of the 
corresponding aircraft, 

e The cumulative time including ground time spent at each stop, without a crew change 
should be less than 24 hours, 

e Every airfield on the route should be compatible with the aircraft, 

e The length of the route should be less than the upper bound* which is defined by the 
user, 

e All airfields can be used by all types of aircraft as an APOD or an APOE but only 
specified aircraft as an enroute, 

e The distance between two adjacent enroute airfields should be greater than a 
prespecified distance defined by the user. 

The most obvious way to find a solution is to devise a bookkeeping scheme that 
allows an orderly representation of all possible paths. It is useful to note that the 
bookkeeping scheme must not allow a node to be revisited (i.e., a cycle) in the network 
With cyclic paths eliminated, networks are equivalent to trees. The tree shown in Figure 
2.3 1s made from the network in Figure 2.2 by using breadth first search and eliminates 
certain routes (e.g, ADBG) by using physical and strategic constraints. Pseudo code for 
the route generation algoritm is described below and additional details are contained in 


Appendix D 





’ Upper bound is defined by a parameter A times the great circle distance between origin and destination. 
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Figure 2.3 A Tree Derived from the Network of Figure 2.2 


L. Pseudo Code of Route Generation Model 


read input data files 
compute the great circle distances between airfields and add a random number between 0-1° 
for dclivery/recovery routes 
for cach aircraft 
detcrmine the adjacency relation between airfields 
if (distance < ferry range ) and (airfield <> airficld2 ) and 
(aircraft is compatible with airfield ) then 
adjacent = 1.0 
otherwise 
adjacent = 0.0 
end if. 
route generation algorithm 
put the origin as first stage in the tree 
LOOP : stage 2 to stage 5 
LOOP : if eof (previous stage) then exit (of LOOP) 
node:= first(previous stage) 
if node # destination then begin 
select nodes which are accessible from node due to 
given constraints 
exclude any node which is found 1n the past route list 
add rest of the nodes to next slage 





* If the geographic location values are not available based on security reason, program can just read the 
distance matrix provided by user. By adding a small random number, the program perturbs the data in 
order to facilitate breaking tics and having routes uniquely identified by their length. 
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end 
goto LOOP 
increment stage level 
end LOOP 


LOOP start from the last stage to first stage 
find the destination node 
follow past route list through origin 
record the route 
end LOOP. 
feasibility checks 
if route length > upper bound then delete the route 
if route is infeasible duc to crew stages then delete the route 
manual! screcning of the routes by user 
cnumerate routes bascd on their length and store for each aircraft 
end for 
manual route adding by user 
end for 
output the route set 
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IV. POSTPROCESSING MODULES 


The final output of the combined centroid and route selection preprocessing step is 
a set of candidate routes that can be used as input data for Throughput II. However the 
airfield aggregation and subsequent route selection process may still generate a large 
number of routes, and it may not be computationally practical to include all the candidate 
routes in the Throughput IT model. Thus, we utilize a method to convert the large 
problem into smaller problems of manageable size. The approach to this reduction is to 
solve a linear program optimized over a subset of its columns, the result of which may be 
suboptimal. Based on ideas from LP sensitivity analysis we can test whether the set of 
decision variables associated with a certain route, not in the original subset, is attractive. If 
it 18 attractive, it can be put into the basis as a new column. Moreover, we can examine 
why this route might be attractive and, if necessary, modify the rules that select the initial 
set of routes in the heuristic route selection model. The Throughput II model may be 


regarded as an example of the following LP model : 


min cx GB 1) 
wh Ama) 
x20 
and its dual 
max 76 (35:2) 


yf 6 WwASe 


mw unrestricted 
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The optimality conditions for a linear program state that necessary and sufficient 


conditions for x° and z° to be an optimal primal-dual pair are, 


1. Primal Feasibility . Ax ~- 6b 
2 Dual Feasibility © 2zA<ec 
3 Strong Duality : 6% =9d 


In standard column generation techniques, new columns are defined with respect 
to the original constraint set. Throughput II, however, does not fall into this framework. 
Instead, associated with each new route-aircraft combination at any point in time is a set 
of decision variables that appear in the original constraints plus an additional set of 
constraints that contain only the new decision variables. For example, one of the new 
constraints specifies that the amount of cargo that can flow along a new route is restricted 
by the capacity of aircraft flying this route. We denote the set of variables associated with 


the new routes y and express their inclusion in model (3.1) via, 


min cx+ fy (3.3) 
st. Ax+hy=b 
Dy =0 
x,y20 


The fact that the right-hand-side of the new set of constraints is zero reflects our 
assumption that the new LP (3.3) is feasible if y=0 ; 1.e., we are not required to fly the 


new routes. The dual of (3.3) is, 
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max 7b (3 4) 
st mAse 
nF +yD< f 


As a result of our non-standard form the usual optimality test regarding the sign of 
the reduced cost of a new column no longer directly applies. However, the following 


theorem provides an analogous optimality test for the framework outlined above. 


Theorem 3.1 


Let x* and 7” be an optimal primal-dual pair for (3.1) and (3.2). If mFi-f,<0 VW, 
then (x, 0) and (z", 0) form an optimal primal-dual pair for (3.3) and (3.4). 
Proof of Theorem 3.1: 
We will verify the Karush - Kuhn - Tucker (KKT) Optimality Conditions of 
1 Primal Feasibility 


2. Strong Duality 
3. Dual Feasibility 


of the primal-dual solution pair (x,y) = (x",0) and (7",y") = (",0) in 


(3.3) and (3.4) 
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1, Primal Feasibility : 


Ax" +hy"=b 

Dy" =0 

x" y">0 
Now, 


Ax’ +l0=6 since Ax ~- 6 and 
DO=0 


(x" 0) = 0 also follows from feasibility of x’ in (3.1). 


2. Strong Duality : 


cx" + fy" = m'b+y"0 
cx’ + f0=n°b+y0 


cx =7°b 


follows from strong duality of x° and 7” in (3.1) and (3.2). 


3. Dual Feasibility : 
n"A,~c, =m A,—e€, <0 follows from dual feasibilty in (3.2). 
Dy = fe =a F, =f, 20 

follows from the optimality test hypothesis. QED 


From Theorem 3.1 we know that if we solve Throughput II and subsequently 


determine that the reduced costs of variables associated with routes not explicitly included 


in the original model pass the optimality test 2° /’, — f, < 0 then the current solution is 
p p y Aes 
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optimal with respect to all routes. If the optimality test is not satisfied then we can append 
the variables associated with attractive routes and re-solve the corresponding Throughput 
Il model. To implement the column generation technique in Throughput I, note that the 
variables associated with the new routes are X, ¥, TONSUE, and 7PAX. Thus, summing 
the dual variables associated with constraints that these variables appear in computes the 
“1°” terms. Additionaly, only the TONSUE and 7PAX variables have objective function 
coefficients (for late penalty computation). Consequently, only these columns have “f” 


terms associated. The complete GAMS code is included in Appendix C. 


It is desirable to prevent, or at least recognize, when Throughput II generates 
solutions that are infeasible with respect to the original, disaggregated network. One 
example of such a violation may occur when a unit’s origin airfield (and hence the unit) is 
aggregated with a centroid closer to the destination airfield; in this case, the origin- 
destination distance may be greater than the selected aircraft’s minimum payload range, 
even though the centroid-destination distance was not longer than the minimum payload 
ferry range. In reality, therefore, this aircraft wold have to use a different route or carry a 
lighter load. Another potential infeasibility can arise when Throughput IJ inappropriately 
“borrow” airfield capacity between airfields that have been aggregated into a single 
airfield. A promising and important area for future research is to implement a 
disaggregation procedure that maps aircraft missions to the original network. With such a 
procedure in place, it will be possible to recognize and take steps to avoid the 


infeasibilities described above. 
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V. RESULTS AND SUMMARY 


This chapter reports the results for model instances based on two different data 
sets provided by AFSAA in the summer and _ fall of 1995. TPFDD data must be 
aggregated for use in Throughput IJ. This aggregation procedure is implemented in the 
Pascal programming language. It takes the input TPFDD, and then performs aggregation 
based ona strategy determined by several user-specified parameters. The flow diagram of 
the aggregation algorithm is in Appendix D and may be summarized as follows: 

Step |: Aggregate units that are to be shipped between the same origin-destination 
pair and have identical delivery time windows based on the available-load-date (ALD) 
and required-delivery-date (RDD) 

Step 2: Eliminate units that have very small movement requirements (e.g., one 
passenger) 

Step 3: Aggregate units that are to be shipped between the same origin-destination 
pair and have similar RDD values. The definition of “similar” is allowed to depend on 
when over the deployment period the units should be delivered. Greater latitude is 
permitted near the end of the scenario compared to the early days of the deployment. This 
is implemented by the user specifying parameters R,, 2, R3, partitioning the planning 
period into segments with corresponding tolerance values r), r2, 73, where, 


Osn<4K<r 


37 


If two units with the same origin and destination have RDD’s in the same segment 
and if the difference between the RDD’s is less than the corresponding r value, the two 


units are aggregated. This is depicted in Figure 5.1 











| RDDI-FRDD2* <P, :RDDI-RDD2! <r, | RDD] -RDD2\< fF, 
L =I — | 
0 Ry Ro R3 








Figure 5 1. Unit Aggregation Over the Time Horizon 


After the aggregation the RDD and ALD values are defined as: 


RDD,..w = Max (RDDI.RDD2) , 

ALD,..w = Weighted Based on Movement Requirements (ALD!,ALD2) 

Step 4: The results of the aggregation procedure are then used in the centroid 
model, followed by the route selection procedure which in turn provides input data for the 


Throughput I] optimization model. 


A. RESULTS 

1. Centroid Model 

The data set that has been used for the centroid models consists of 117 airfields in 
the Continental United States (CONUS). All the models presented in Chapter III can be 
used to select the centroid airfields. In order to investigate the relative merits of these 
models three criteria have been defined. The first one is the computational time required to 


solve the model on an IBM RS6000/590 system using GAMS with the OSL solver. The 
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second criterion is the average distance between the centroid airfields and their satellite 
airfields. This is an important issue in order to have an accurate delivery cycle time in the 
Throughput IJ optimization model. The last criterion is the primary objective value of the 
proposed model which 1s the maximum distance between the centroid airfields and their 
satellite airfields. The model that we develop allows the user to predefine certain centroid- 
satellite groups and effectively remove them from consideration in the centroid model. 
This ts desirable feature because it allows the user to prevent the possibility of capacity 
borrowing at known system bottlenecks. The results are illustrated in Table 5.1. 


TIME AVERAGE 
(RS6000/590 DISTANCE 


MAXIMUM 
DISTANCE 


GAMS/OSL) (nm) (nm) 


YOST’S MODEL 


14 SECONDS 





MINISUM MODEL 


10 MINUTES 








MINIMAX 
MODEL* 


32 SECONDS 





SUGGESTED 
MODEL 


Table 5.1. Computational Results for Various Centroid Models. 


87 SECONDS 














Note that Yost’s model is a minisum model with a limited number of airfields that 


may be used as centroids. As a result, in runs quickly but performs worse than the other 


models relative to the other criterion. The suggested model provides an attractive balance 


between the minisum model and the minimax model with respect to the average and 


maximum distance criterion and requires modest computational effort. 





° This is the modcl that solved via sct covering algorithm. 


Recail that the centroid models require the user to define the (maximum) number 
of centroids that may be selected. Figure 5.2 illustrates how the maximum centroid- 
satellite distance varies with the number of centroids. As shown in Figure 5.2, there is a 
large distance when the centroid budget is small. Large range values are more likely to 


lead to infeasible (unflyable) routes in the route selection model. 
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Figure 5.2. Computational Results for Centroid Model. 
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z Route Selection Model 


After we select the centroid airfields for origin and destination airfields via the 
centroid model, the route generation model constructs routes based on user-defined input 
parameters. It takes less than an hour to set up the input files for the first run. After this 
preparation, each run takes approximately five minutes. 

The first input parameter is denoted 4 and is used to define the maximum route 
length between origin-destination pairs. This lambda value actually consists of the product 
of three different parameters. We assign a lambda value for each airfield and aircraft. To 
illustrate this process, if we call the resultant lambda value as the composite lambda value, 
then, 


* A 


“composite — ““ ongin destination 


A 


aircraft 


If we define each of these parameters as 1.0, then only the direct route (i.e., no 
enroute stops) will be generated. Using these three different parameters provides flexibility 
for selecting routes based on the type of aircraft or accessibility of different regions of the 
world with origin/destination dependent lambda values 

The second parameter of the route selection model is the restriction on the flight 
time of an aircrew. In this research we used 24 hours for the maximum flight time without 
a crew-change. The final input required from the user is the set of crew-stage airfields and 


enroute airfields that can be used during the airlift operation. Table 5.2 shows the number 


4] 


of infeasible routes for various A values for the C747 which has a minimum payload 


ferry range of 6500 nm for delivery routes and 9200 nm for recovery routes, 







Total Routes 





Composite 


Lambda Value 


Total 


Feasible Routes 





Total Infeasible Routes Due 


To Flight Time Restrictions 





| 





7266 





1.041 


LULZ 











Delivery Routes 1.161 
1.184 

1.240 

1.000 

4152 1.139 
Recovery Routes 1.162 


1.185 





233, 





Table 5.2. Route Selection Model Results. 
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We have performed experiments in order to show the “quality” of the routes that 


have been generated by the route selection procedure versus the routes with the same data 


set in Lim (1994). The computational results based on comparing manually generated 


routes and the model are shown in Table 5.3. 
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f T 
Model Number of Objective Function 
Routes | Value 
Route Selection Model 71 routes 13.208 
Manually Generated 62 routes 15.021 
Routes 











Table 5 3. Route Selection Model’s Results. 


The route selection model improves the routing options in Throughput II so that 


there is a 12% reduction in the weighted sum of late and undelivered cargo. 


3. Postprocessing Module 


One important issue in the postprocessing module is whether including additional 
routes in Throughput II will result in an improved solution. There are many decision 
variables associated with a single physical route through the network: 

The route may be flown by several types of planes and at different points in time in 
the planning period. There are typically a large number of decision variables associated 
with routes not explicitly included in Throughput II that have attractive reduced costs (see 
Chapter II). However, many of these provide small if any improvement in the solution, if 
included. For this reason, we only augment the set of routing options with decision 
variables whose reduced costs are significantly attractive, relative to some tolerance. In 
the proposed method, Throughput J] is solved with respect to the original routes and then 


the set of routes is augmented and the model is re-optimized. 
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When this is done, it is desirable to initialize the re-optimization with the basis 
from optimal solution of the original model. Table 5.4 illustrates the savings in 
computational effort that this technique provides relative to solving the larger model from 
scratch. In addition, Table 5.4 shows that the postprocessing scheme found routing 


options that significantly decreased the weighted sum of late and undelivered cargo. 











Number of Model Objective | Running Time 
Routes Function Value ___(RS6000) 
71 routes Throughput II 13.208 10 — minutes 








8.0292 60 minutes 


8.0292 38.75 minutes 


Table 5.4. Computational Results of Sensitivity Analysis 


117 routes Throughput IT 












117 routes Proposed Model 











B. RECOMMENDATIONS AND FUTURE WORK 


The techniques proposed in this thesis have been recently applied to a large 
TPFDD provided by AFSAA in August 1995 for a C17/C747 fleet mix study. Due to the 
enormous size of the data, which has approximately 14500 movement requirements, the 
preprocessing modules for centroid selection and route reduction were an invaluable tool 
for quickly generating an optimization model of tractable size. This modeling and data 
reduction exercise also revealed that certain additional features could further streamline 
the process for future analyses: 

e User interface: Current aggregation and route selection modules are implemented as a 


Disk Operating System (DOS) application. Because of the structure of this operating 
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system, the ease of use for an analyst would be greatly enhanced by implementing this 
software as a Windows application with user-friendly interfaces. 

e Data Interface: The current model uses text files that must be manually entered by the 
analyst. Because this is a time consuming procedure, a data interface that permits easy 
modifications while ensuring data integrity would significantly decrease the model 
generation time. 

e Route Generation Module: As indicated above, there is a set of decision variables 
associated with each route. It would be useful to develop a theory that utilizes reduced 
cost information from individual decision variables to indicate when the inclusion the 
set of decision variables will improve (or be likely to improve) the current solution. 

e Feasibility Test Postprocessing Module: As discussed at the end of Chapter IV, 
deriving a procedure that maps Throughput II solutions to the original disaggregated 
network would prove valuable with respect to testing the feasibility of the proposed 
solution 

As shown above, much work remains to realize the full potential of pre and post- 


processing of the Throughput IT model. The contributions of this thesis form a solid 


foundation for these improvments. 
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APPENDIX A. THROUGHPUT I MODEL 


The following is a brief summary of Throughput II (Morton e¢. al, 1995) 


A. 


INDICES 

u indexes units, e.g., 82nd Airborne 

a indexes aircraft types, e.g., C5, C141 

ie index time periods 

b indexes all airfields (origins, enroutes and destinations) 
i indexes origin airfields 

k indexes destination airfields 

r indexes routes 

INDEX SETS 

1, Airfield Index Sets 

B set of available airfields 

[cB origin airfields 

KcR destination airfields 

23 Aircraft Index Sets 

A set of available aircraft types 

A buk&A aircraft capable of hauling bulk size cargo 
AoverCAbu aircraft capable of hauling over-sized cargo 
AwwCAwe — aircraft capable of hauling out-sized cargo 


Bulk cargo is palletized on 88 x 108 inch platforms, which can fit on any plane. 


Over-sized cargo is non-palletized rolling stock, it is larger than bulk cargo and can fit on 


a C141, C5 or C17. Out-sized cargo is very large non-palletized cargo that can fit into a 


C5 or C17 but not a C141. 
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3. Route Index Sets 


R set of available routes 

RacR permissible routes for aircraft a 

RupcRa permissible routes for aircraft a that use airfield 5 

RarcRa permissible routes for aircraft a that have origin # and destination k 
DR. CR delivery routes that originate from origin / 

RRR recovery routes that originate from destination k 

4. Time Index Set 

if set of time periods 

Ye possible launch times of sorties for unit « using aircraft a and router 


The set 7),., covers the allowed time window for unit w, which starts on the unit’s 
available-to-load date and ends on the unit’s required delivery date, plus some extra time 


up to the maximum allowed lateness for the unit. 


Cc. GIVEN DATA 


1. Movement Requirements Data 


MovePAX ux Troop movement requirement for unit u from origin i to 
destination k 


Move UF Equipment movement requirement in short tons (stons) for unit w 
from origin / to destination k 


Probulk,, Proportion of unit w cargo that is bulk-sized 

ProOver, Proportion of unit w cargo that is over-sized 

ProOut, Proportion of unit v cargo that is out-sized 

2. Penalty Data 

LatePenUE,, Lateness penalty (per ston per day) for unit « equipment 
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LatePenPAX, Lateness penalty (per soldier per day) for unit w troops 
NoGoPenUE,, Non-delivery penalty (per ston) for unit « equipment 
NoGoPenPAX, Non-delivery penalty (per soldier) for unit u troops 
MaxLate Maximum allowed lateness (in days) for delivery 


Preserve a Penalty (small artificial cost) for keeping aircraft a in mobility 
system at time / 


3: Cargo Data 

UESQI 1, Average cargo floor space (in sq. ft.) per ston of unit « equipment 

PAXW1, Average weight of a unit w soldier inclusive of personal 
equipment 

4. Aircraft Data 

Suppbya Number of aircraft of type a that become available at time / 

MaxPAX, Maximum troop carriage capacity of aircraft a 

PAXSGF tua Average cargo space (in sq. ft.) consumed by a unit w soldier fro 
aircraft a 

ACSq!1, Cargo floor space (in sq. ft.) of aircraft a 

Loadl ff, Cargo space loading efficiency (<1) for aircraft a. This accounts 


for the fact that it is not possible in practice to fully utilize the 
cargo space. 


URate, Established utilization rate (flying hours per aircraft per day) for 
aircraft a 

S. Airfield Data 

MOGCap, Aircraft capacity (in narrow-body equivalents) at airfield 6 
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MOGRequs 


Conversion factor to narrow-body equivalents for one aircraft of 
type a at airfield 5 


MOGISf MOG efficiency factor (<1), to account for the fact that it ts 
impossible to fully utilize available MOG capacity due to 
randomness of ground times 

6. Aircraft Route Performance Data 

Maxload,, Maximum payload (in stons) for aircraft a flying route r 

Gime ahr Aircraft ground time (due to onload or offload of cargo, 
refueling, maintenance, etc.) needed for aircraft a@ at airfield 6 on 
route r 

DTimean, Cumulative time (flight time plus ground time) taken by aircraft a 
to reach airfield 4 along route r 

PtTimear Total flying hours consumed by aircraft a on route r 

CTimeo, Cumulative time (flight time plus ground time) taken by aircraft a 


DaysLate yar 


on route r 


Number of days late unit u’s requirement would be if delivered by 
aircraft a via route r with mission start time / 


DECISION VARIABLES 


1. Sortie Variables 


Xan 


Yan 


Number of aircraft a that airlift unit w via route r with mission 
start time during period ¢ 


Number of aircraft a that recover from a destination airfield via 
route r with start time during period ¢ 
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2. Aircraft Allocation and De-allocation Variables 


Allotar Number of aircraft a that become available at time / that are 
allocated to origin / 


Releases Number of aircraft a available at origin / in time / that are not 
scheduled for any flights from then on 


a Aircraft Inventory Variables 

Aan Number of aircraft a inventoried at origin / at time / 

AP aks Number of aircraft a inventoried at destination & at time / 
NPlaneSar Number of aircraft a in the air mobility system at time / 

4. Airlift Quantity Variables 

Tons UE yan Total stons of unit w equipment airlifted by aircraft a via route r 


with mission start time during period / 


TPAX vant Total number of unit u troops airlifted by aircraft a via route r 
with mission start time during period / 


5. Elastic (Nondelivery) Variables 


UL NoGouk Total stons of unit w~ equipment with origin / and destination k 
that is not delivered in the prescribed time frame 


PAXNoGo,«. Number of unit uw troops with origin / and destination k who are 
not delivered in the prescribed time frame 
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E. OBJECTIVE 
minimize 


SDD. Dd LatePenvl:, * DaysLate,,,,* TonsuE (A.1) 


ua reR, teTay, 


+> YY LatePenPAX ,*DaysLate,,,,* TPAX 


uo a ER Tugy 


+) ¥ (NoGoPenvE,,* VENoGo,, + NoGoPenPAX , * PAXNoGo,, ) 
u 1 k 


uart 


uart 


+)° > Preserve,,* NPlanes,, 
a t 


The objective function minimizes the total weighted penalties incurred for late 
deliveries and non-deliveries. The model’s secondary objective is to choose a feasible 


solution that maximizes unused aircraft. 


F. CONSTRAINTS 


> dX DY Tonsve,,,, + VENoGo,, = MoveVE,,, (A.2) 
a Angie 7 Regt 1 Toy 
Vu,i,k: MoveUE,, > 0 
> & & Tonsvr.,,, + VENoGo,, 2 ProOut,* MoveUE,,, (A.3) 
a Aan TERae OT 
Vu,i,k: MoveVE,, > 0 
» SY VonsUE 4, + UENOGO,,, 2 (A.4) 
UC Ag, TORigg Tay . 
(ProOver, + ProQut,|* MoveVE,,, u,i,k: MoveVE,,, > 0 
YD Dd TPAX + PaxNoGo,, =MovePAXx (A.5) 


a TOR ae ECT ay 


V u,i,k: MovePAX,, >0 
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X +H. + Release, = (A 6) 
> > uart aul aut 


r oreDR, 


H,, , + Allot,,+ > Sle x WGA 


re Ry CspCTimeg} 
> ¥, + HP, = HF y , 1 +>) > be ae Va,k,t (A.7) 
rv RR, u Nea sitet det 
> >, Allot,,, < Y supply, Val (A.8) 
role gry 
NPlanes,, = >- > Allot,,. — >. > Release,,., Val (A.9) 


wou? a i) 


Zs ee ¥ Bg + 2, > Kg? Yor - ee Hay (A.10) 


OR UL oa rR, UN ear 


2" HP. < y NPlanes,,, Val 
4 


ae oe | 
where 
r (-t' +1, if t’st<t'+CTime, -1 
“| CTime,,, if 12> t'+CTime,, -1 


TPAX ,,,, < MaxPAX ,* X 4, VujartuteT,, (A.11) 

TonsUL: ,, + PAXWt* TPAX |, < MaxLoad,,* X ,,,, (A.12) 
Vugar.tteT,, 

PAXSqIt,.* TPAX ,,,, + VESGFt, * TonsuE,,, S (A.13) 


ACSqlt, * Loadlff* X Via treT, 


uart > 
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> YD Mette, * X jan + >> Fuvime,, * | es 
reR, t 


eR Tay 


> URate,* NPlanes,,, Va 


YS DY (MOCRe|,,* GTime,,, 1 24)* X47 


ua rR, "eT 


me Timi, | t 

+e > (MOGRegq,,*GTime,,, / 24)*¥,,, 
aorcR, t'+[ DTimeg, |=! 

< MOGEA * MOGCap, , Vout 


A.2 Demand satisfaction constraints for all classes of cargo 
A.3 Demand satisfaction constraints for out-sized cargo 
A.4 Demand satisfaction constraint for over-sized cargo 
A.5 Demand satisfaction for troops 

A.© Aircraft balance constraints at origin airfields 

A.7 Aircraft balance constraints at destination airfields 


A.8 Aircraft balance constraints for allocations to origins 


A.9 Aircraft balance constraints accounting for allocations and releases 


A.10 Cumulative aircraft balance constraints 
A.|1 Troop carriage capacity constraints 
A.12. Maximum payload constraints 

A.13 Cargo floor space constraints 

A.14 Aircraft utilization constraints 


A.15. Airfield MOG constraints 
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(A.14) 


(A.15) 


APPENDIX B. GAMS CODES FOR CENTROID MODEL 


$TITLE SETCOVERING 
$STITLE Yasin TURKER 29 JUL 1995 


SOFFUPPER OFFSYMLIST OFFSYMXREF INLINECOM { } 
$ONMIXED 
OPTIONS 
LIMCOL = 0 , LIMROW = 0). , SOLPRINT = OFF , DECIMALS = 2 
RESLIM = 18000, ITERLIM = 400000, OPTCR =0.1,SEED =3141 
MIP = OSL 
*include centroid model for origins 
$INCLUDE origset.gms 
SETS INDEX /1*100/ 
AF airficlds 
$INCLUDE airfield.sct 
\(AF) origins 
$INCLUDE origin] .txt 


ALIAS (1.1P.J): 
ALIAS(AF,AFP.AFPP); 


PARAMETER DIS (AF.AFP) 
$INCLUDE dist] .txt 


SCALAR K2_ number of airficlds which is allowed throughput model /4/ 
RAD | radius of centroid 
RBEST best radius of centroid 
UPP upper value of centroid 
LOW lower value of centroid; 
UPP=SMAX((1,J),DIS(I,J)), 
LOW=SMIN((1,J).DIS(.J)), 
RAD=(UPP+LOW)/2, 
RBEST=UPP, 


VARIABLE 
ZORIG maximum number of centroids 
ZIORIG objective function value of minisum model, 


BINARY VARIABLE 
YORIG(J) open centroid j; 


EQUATIONS 
OBJ objective function of covering problem 
OB]JI objective function of minisum model 


COVER(|) covering assignments for sct covering model 
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ASSIGNI(I.IP) assignment for cach airfield for an centroid 
ASSIGN2 number of airficld which is allowed for throughput model 
ASSIGN3(1) make sure every airficld either be in centoid or alone, 


*MINIMIZE 
OBJ.. SUM(J,YORIG(J))=E=ZORIG,; 
COVER(1).. SUM, YORIG(J)$(DIS(1,J) LE RAD))=G=1,; 


MODEL ORIGCOV /OB],COVER’, 


LOOP(INDEX$(NOT(((UPP - LOW) LE 5.0) AND (ZORIG.L EQ K2))), 
RAD=(UPP+LOW)/2; 
SOLVE ORIGCOV USING MIP MINIMIZING ZORIG, 
DISPLAY RAD.UPP.LOW; 
UPP=RADS$(ZORIG.L LE K2+0.1)+UPP$(ZORIG.L GE K2+0.9); 
LOW=RADS$(ZORIG.L GE K2+0.9)+LOW$(ZORIG.L LE K2+0.1); 


): 
DISPLAY YORIG.L.RAD,UPP,LOW,ZORIG.L, 


POSITIVE VARIABLE 
XORIG(1.IP) airficld i belongs to centroid j ; 


XORIG.UP(I,J)$(DIS(.J) LE RAD) =1.0; 


*minisum modcl 

OBJ1.. SUM((1,J))$(DIS(1.J) LE RAD),DIS(1,J)*XORIG(I,J))=E=Z1ORIG; 
ASSIGN 1(1,J)$(DIS(I,J) LE RAD).. XORIG(I,J)=L=YORIG(J); 

ASSIGN2.. SUM(J, YORIG(J))=E=K2; 

ASSIGN3(1).. SUM(J$(DIS(I,J) LE RAD), XORIG(I,J)) =E=1,; 

MODEL CENTROID/OBJ1,ASSIGNI,ASSIGN2,ASSIGN3/, 

SOLVE CENTROID USING MIP MINIMIZING Z1ORIG; 

DISPLAY YORIG.L,ZORIG.L,ZIORIG.L,XORIG.L; 


*include centroid modcl for destinatinations 
$INCLUDE destsct.gms 

SETS K(AF) airports 

$INCLUDE destinal.txt 


ALIAS (K,KP.L); 
UPP=SMAX((K,L),DIS(K.L)); 
LOW=SMIN((K,L),DIS(K.L)): 
RAD=(UPP+LOW)/2: 


SCALAR K1 number of centroid destination airfields is allowed in Throughput /8/; 
VARIABLE 

ZDEST — maximum number of centroids 

ZIDEST objective function value of minisum model, 


BINARY VARIABLE 
YDEST(L) open centroid L, 
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EQUATIONS 

OBJ4 — objective function of covering problem 

OBJ6 — objective function of minisum model 

COVERI(K) covering assignments for set covering model 
ASSIGNII(K.KP) assignment for cach airfield for an centroid 
ASSIGN21 number of airficld which is allowed for throughput model 
ASSIGN31(K) make sure every airfield cither be in centoid or alone; 


*MINIMIZE 
OBJ4. SUM(L,YDEST(L))=E=ZDEST, 
COVERI(K).. SUM(L,YDEST(L)$(DIS(K.L) LE RAD))=G=1; 


MODEL DESTCOV /OBJ4.COVERI/; 


LOOP(INDEX$(NOT(((UPP - LOW) LE 5.0) AND (ZDEST.L EQ K1))), 
RAD=(UPP+LOW)/2, 
SOLVE DESTCOV USING MIP MINIMIZING ZDEST; 
UPP=RAD$(ZDEST.L LE K]+0.1)+UPP$(ZDEST.L GT K1+0.9); 
LOW=RADS$(ZDEST.L GT K1+0.9)+LOWS(ZDEST.L LE K1+0.1); 


yd 
DISPLAY YDEST.L,RAD.UPP,LOW.ZDEST.L,; 


POSITIVE VARIABLE 
XDEST(K,KP) airfield i belongs to centroid L ; 
XDEST UP(K.L)$(DIS(K,L) LE RAD) =1.0, 


*minisum modcl 

OBJ6.. SUM((K.L) $(DIS(K,L) LE RAD),DIS(K,L)*XDEST(K,L))=E=Z1 DEST, 
ASSIGNI1(K,L)$(DIS(K,L) LE RAD).. XDEST(K,L)=L=YDEST(L), 
ASSIGN21.. SUM(L,YDEST(L))=E=K 1, 

ASSIGN3I(K).. SUM(L$(DIS(K,L) LE RAD), XDEST(K,L)) =E=1, 

MODEL CENTROID1I/OBJ6,ASSIGN11,ASSIGN21,ASSIGN31/; 

SOLVE CENTROIDI USING MIP MINIMIZING Z1DEST, 

DISPLAY YDEST.L,ZDEST.L,Z1DEST.L,XDEST.L; 


*include sct of units 

SET UNIT units in TPFDD 
/ 

$INCLUDE unit] txt 

Es 


SETS AB aircraft body type 
*Auircraft arc classified as WB (widcbody), NB (narrow body), 
*GM (ground manuverablc), or TT (tactical). 

/ WB, NB, GM, TT / 


ENR(AF) cnroute airfields 


$INCLUDE enroute.txt 
i 
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*moving cargo requirements for units given origin and destination 
*origins and destinations arc not centroids 

PARAMETER TONSUE(UNIT,I.K) 

$INCLUDE moving ].txt 


*moving passcnger requirements for given origin and destination 
*origins and destinations arc not centroids 

PARAMETER TPAX(UNIT,I,K) 

$INCLUDE movepa! 1.txt 


*mog capacity for all airficld before aggregation 
TABLE MOG(AF,AB) 
$INCLUDE mogcap.txt 


scalar totton total cargo aficr aggregation 
tottpax total passenger after agercgation; 

totton=sum((unit,i,k),tonsue(unit,i,k)); 

tottpax=sum((unit,1,k),tpax(unit,1,k)); 


*apercgation parameter 

*ageregatc cargo and passenger duc to its origin and centroid 

PARAMETER REPORT2(UNIT,J,L), 

PARAMETER REPORT22(J,L); 

REPORT 2(UNIT,J.L)=SUM((I,K)$((XORIG.L(I,J) EQ 1.0) AND (XDEST.L(K,L) EQ 1.0)), 
TONSUE(UNIT,I,K)); 

REPORT22(J,L)=SUM(UNIT,REPORT2(UNIT,J.L)); 


PARAMETER REPORT3(UNIT,J,L); 

REPORT3(UNIT,J,L)=SUM((I,K)$((XORIG.L(I,J) EQ 1.0) AND (XDEST.L(K,L) EQ 1.0)), 
TPAX(UNIT,I,K)); 

PARAMETER REPORT6(J,AB), 

PARAMETER REPORTG61(L,AB), 

*apercgation of mog valucs into centroids first origin 

REPORT6(J.AB)=SUM(I$(XORIG.L(1,J) EQ 1.0). MOG(1,AB)); 

*now destinations 

REPORT61(L,AB)=SUM(K$(XDEST.L(K,L) EQ 1.0),MOG(K,AB)); 


PARAMETER REPORTA(J.L), 
LOOP(UNIT, 

REPORT4(J,L)$((REPORT2(UNIT,J,L) GT 0.0) OR (REPORT3(UNIT,J,L) GT 0.0))=1, 
); 
PARAMETER REPORTS(*); 
REPORTS(J)$((SUM(L,REPORT4(J,L))) GT 0.0 )=1.0; 
REPORT S(L)$((SUM(J.REPORT4(J,L))) GT 0.0) =1.0; 


FILE Fl /movement.txv/; 
FILE F2 /distance.txt/; 
FILE F3 /pairs.txv/: 
FILE F4 /mogcap.dav; 
FILE FS /origin1] 1.txv/; 
FILE F6 /destin! 1.txt/, 
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FILE F7 /distance.dav/. 
FILE F8 /routclog.txv. 


PUT FI. 
OPTION REPORT2:2:0:1. 
OPTION REPORT3:2:0:1; 
PUT "MOVEMENT REQUIREMENTS". 
LOOP((UNIT,J.L)$((REPORT2(UNIT,J.L) GT 0.0) OR REPORT3(UNIT,J,L) GT 0.0), 
PUT 7: 
PUT UNIT.TL:4. 
PUT J TL:4; 
PUT LoThL4: 
PUT* 
PUT REPORT2(UNIT,J.L), 
PU ot 
PUT REPORT3(UNIT,J,L); 
yi 


PUT F2. 
PUT "DISTANCES BETWEEN AIRFIELDS", 
LOOP ((AF,AFP)$((ENR(AF) OR I(AF) OR K(AF)) AND 
(ENR(AFP) OR I(AFP) OR K(AFP))) , 
PUT /. 
PUT AF.TL:4, 
PUT AFP. TL:4, 
PUT” -*: 
PUT DIS(AF.AFP). 
i 


PUT. F3? 
PUT “origin destinauon pairs for route selection model"; 
LOOP((J.L)$(REPORT4(J,L) EQ 1.0), 

PUT /, 

PUT J.TL:4; 

PUT L.TL:4,; 


). 


PUT F4. 
PUT " WB NB GM TT"; 
LOOP(ENRS$(REPORTS(ENR) NE 1.0), 
PUT /, 
PUT ENR.TL:4, 
PUT” 3 
PUT MOG(ENR,'WB)):2:0; 
PUT Ns 
PUT MOG(ENR.'NB’):2:0; 
PUT” “¢ 
PUT MOG(ENR,'GM!):2:0; 
PUT" "3 
PUT MOG(ENR,'TT’):2:0: 





LOOP(J$(REPORTS(J) EQ 1.0), 
PUT /. 
PUT J.TL:4, 
PUR 
PUT REPORT6(J,'WB'):2:0, 
PUES *s 
PUT REPORT6(J.'NB'):2:0; 
Pus 3 
PUT REPORT6(J,'GM'):2:0; 
Pur 
PUT REPORTO6(J.'TT'):2:0: 
). 





LOOP(L$(REPORT5(L) EQ 1.0), 
PUT /. 
PUT L.TL4: 
PUT 
PUT REPORTGI(L,'WB'):2:0, 
Puy” 
PUT REPORTGI(L,'NB’):2:0: 
POLE #5 
PUT REPORT6I(L,'GM'):2:0; 
eur” 
PUT REPORT6I(L,'TT’):2:0, 
) 





PUT FS; 
PUT "ORIGIN CENTROID AIRFIELDS", 
LOOP(J$(YORIG.L(J) EQ 1.0), 
PUT /; 
PUT J.TL:4, 
R 


PUT F6, 
PUT "DESTINATION AIRFIELDS", 
LOOP(L$(Y DEST.L(L) EQ 1.0), 
PUT /; 
PUT L.TL:4, 
) 


PUT F7, 
PUT "*DISTANCES BETWEEN AIRFIELDS"; 
PUT /. 
PUT "/": 
LOOP ((AF.AFP)8((ENR(AF) OR I(AF) OR K(AF)) AND 
(ENR(AFP) OR I(AFP) OR K(AFP))) , 
PUT /, 
PUT AF.TL:4; 
i alee 
PUT AFP.TL.4; 
purs'= * 
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PUT DIS(AF.AFP). 
). 


PUT /. 
BUT 2". 


PUT F8, 

F8.AP=1, 

PUT "CENTROID MODEL RESULTS", 
PUT /: 

PUT “DATE” “: 

PUT SYSTEM.DATE, 

PUT /: 

PUT "TIME". 

PUT SYSTEM.TIME,; 

PUTS 

PUT "ORIGIN CENTROID AIRFIELDS", 
PUT /, 





LOOP(J$(YORIG.L(J) EQ 1.0), 
PUT TL. TLiae 
PUT /. 

% 


PUT "DESTINATION CENTROID AIRFIELDS"; 
PUT /. 


LOOP(L$(YDEST.L(L) EQ 1.0), 
PUT L.TL:4; 
PUT /: 

% 


PUT "MOG VALUES FOR THE AIRFIELDS", 
PUT /: 
PUT" WB NB GM TT", 
LOOP(ENRS(REPORTS(ENR) NE 1.0), 

PUT /, 

PUT ENR.TL:4, 

port * 

PUT MOG(ENR,'WB’):2:0 

PUT” ", 

PUT MOG(ENR,'NB'):2:0, 

PUT" "; 

PUT MOG(ENR,'GM’):2:0; 

BuT™ * 

PUT MOG(ENR,'TT’):2:0; 

hs 


LOOP(J$(REPORTS(J) EQ 1.0). 
PUT /. 





PUT Jobat: 





PUP 3 
PUT REPORT6(J,"WB'):2:0; 
PUY 
PUT REPORTO(J.'NB'):2:0; 
PUPS 
PUT REPORTO(J,'GM’):2:0; 
POE Ms 
PUT REPORT6(J.'TT'):2:0: 
). 


LOOP(L$(REPORT 5(L) EQ 1.0). 
PUT /, 
PUT L.TL:4, 
PUT" 
PUT REPORT61(L.'WB'):2:0; 
PUY % 
PUT REPORT61(L,'NB’):2:0; 
Pur" ". 
PUT REPORT61(L.'GM’):2:0. 
PUT 
PUT REPORT6I(L,'TT'):2:0; 

): 





DISPLAY ORIG.L,XDEST.L,TONSUE,TPAX,REPORT2,REPORT3,REPORT4, REPORT6,REPORT61, 
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APPENDIX C. SENSITIVITY ANALYSIS FOR ROUTES 


$TITLE COLUMN GENERATION 
$STITLE Yasin TURKER 29 APR 1995 


SOFFUPPER OFFSYMLIST OFFSYMXREF 


OPTIONS 
LIMCOL = 0, LIMROW = 0. , SOLPRINT = OFF , DECIMALS = 2 
RESLIM = 18000, [TERLIM = 400000, OPTCR =0.1,SEED =3141 


*read candidate routes 
SETS 
$INCLUDE gamscol.dat 


DELR(GEN) additional routes for column generation 
DDR(GEN) 

DRR(GEN) 

DRI(GEN.AF) 

DRK(GEN,AF) 

DRE(GEN,AF) 

DADJACENT(GEN.AF,AF) 

DRA(GEN.A): 


DELR(GEN)$(CAN(GEN))=YES$(NOT R(GEN)); 

DDR(GEN)$(DRCAN(GEN))=YES$(NOT DR(GEN)); 

DRR(GEN)$(RRCAN(GEN))=YES$(NOT RR(GEN)); 
DRI(GEN,AF)$(RICAN(GEN,AF))=YES$(NOT RI(GEN,AF)); 

DRK(GEN.AF)$(RKCAN(GEN, AF))=YES$(NOT RK(GEN,AF)); 
DRE(GEN,AF)$(RECAN(GEN,AF))=YES$(NOT RE(GEN,AF)); 
DADJACENT(GEN, AF, AFP)$(ADJCAN(GEN, AF, AFP))=YES$(NOT ADJACENT(GEN, AF,AFP)); 
DRA(GEN,A)$(RACAN(GEN,A))=YES$(NOT RA(GEN,A)); 


MAXLEG(GEN)$DELR(GEN)=SMAX((AF, AFP)$DADJACENT(GEN, AF,AFP),DIST(AF,AFP)); 


LOOP((GEN,A)$(DELR(GEN) * DRA(GEN,A)), 
CURRENT(AF)=NO,; 


* current is sct to the first airficld on the route 
CURRENT(AF)=YES$((DDR(GEN) * DRI(GEN,AF)) + (DRR(GEN) * DRK(GEN,AF)) ): 
NEXT(AF)=NO, 
NEXT(AF)=SUM(CURRENT.DADJACENT(GEN.CURRENT,AF)); 
FLTTIME(A.GEN)=SUM((CURRENT,NEXT),DIST(CURRENT,NEXT))/(100*SPD(A)), 
WI(AF)=CURRENT(AF)*I(AF); 
WK(AF)=CURRENT(AF)*K(AF), 
CTIME(A.GEN,NEXT)=GTIME(A,'ONLD')$(CARD(WI) GT 0) + 

GTIME(A,'OFFLD')$(CARD(WK) GT 0) + 
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SUM(CURRENT,DIST(CURRENT,NEXT))/SPD(A); 
KEND=1, 
LOOP(AFP$KEND, 
CURRENT(AF)=NO, 
CURRENT(AF)=Y ES$NEXT(AF): 
NEXT(AF)=NO, 
NEXT(AF)=SUM(CURRENT, DADJACENT(GEN,CURRENT,AF)); 
FLTTIME(A,GEN)=FLTTIME(A,GEN)+ 
SUM((CURRENT,NEXT),DIST(CURRENT,NEXT))/(100*SPD(A)), 
CTIME(A,GEN,NEXT)=SUM(CURRENT,CTIME(A,GEN,CURRENT))+ 
GTIME(A.'ENR')+ 
SUM(CURRENT,DIST(CURRENT,NEXT))/SPD(A); 
* quit inncr loop if NEXT is empty 
KEND= CARD(NEXT), 
) 
): 


MAXLOAD(A.GEN) $ (DELR(GEN) * DDR(GEN) * DRA(GEN,A)) = 
ACLOAD(A,'l')$(MAXLEG(GEN) LE ACRANGE(A,'1')) + 
(ACLOAD(A.'1')+SLOPE(A,'2')*(ACRANGE(A,'I')-MAXLEG(GEN))) 
$(MAXLEG(GEN) GE ACRANGE(A,'I') AND 
MAXLEG(GEN) LE ACRANGE(A,'2')) + 
(ACLOAD(A.'2')+SLOPE(A,'3')*(ACRANGE(A,'2')-MAXLEG(GEN))) 
$(MAXLEG(GEN) GE ACRANGE(A,'2') AND 
MAXLEG(GEN) LE ACRANGE(A,'3')) + 
(ACLOAD(A,'3')+SLOPE(A, '4')*(ACRANGE(A,'3')-MAXLEG(GEN))) 
$(MAXLEG(GEN) GE ACRANGE(A,'3') AND 
MAXLEG(GEN) LE ACRANGE(A,'4')) + 
(ACLOAD(A,'4")+SLOPE(A,'S')*(ACRANGE(A,'4")-MAXLEG(GEN))) 
$(MAXLEG(GEN) GE ACRANGE(A,'4') AND 
MAXLEG(GEN) LE ACRANGE(A,'S')) + 
(ACLOAD(A.'S')+SLOPE(A,'6')*(ACRANGE(A,'S')-MAXLEG(GEN))) 
$(MAXLEG(GEN) GE ACRANGE(A,'S') AND 
MAXLEG(GEN) LE ACRANGE(A,'6')) + 
(ACLOAD(A,'G')+SLOPE(A.'7')*(ACRANGE(A,'6')- MAXLEG(GEN))) 
$(MAXLEG(GEN) GE ACRANGE(A,'6') AND 
MAXLEG(GEN) LE ACRANGE(A,'7’)) 
$ONMIXED 


* first delivery routes 
KCONSUME(A,GEN,TP.T)$( (ORD(TP) LE ORD(T)) AND DELR(GEN) AND 
(DDR(GEN)*DRA(GEN,A)) )=SUM(K$DRK(GEN,K), CTIME(A,GEN,K)/24); 


KCONSUME(A,GEN.TP,T)$( (ORD(TP) LE ORD(T)) AND DELR(GEN) AND 
(DDR(GEN)*DRA(GEN,A)) )=MIN( ORD(T)-ORD(TP)+1 , KCONSUME(A,GEN,TP,T)); 


* now recovcry routcs 


KCONSUME(A,GEN,TP,T)$( (ORD(TP) LE ORD(T)) AND DELR(GEN) AND 
(DRR(GEN)*DRA(GEN,A)) )=SUM(I$DRI(GEN,I), CTIME(A,GEN,1)/24); 
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KCONSUME(A.GEN.TP.T)$( (ORD(TP) LE ORD(T)) AND DELR(GEN) AND 
(DRR(GEN)*DRA(GEN.A)) ) =MIN( ORD(T)-ORD(TP)+I . KCONSUME(A,GEN,TP,T)), 


Hepa tte ttt +++4++4+4++ DYNAMIC SET DECLARATION +++4+++++++H-44+4++4+++ 


DSETX(U,A.GEN,T) is used to control the allowable combination 

of unit. aircraft. route, and start time for the following 

decision variables: X(U,A,GEN,T), TONSUE(U,A,GEN,T) and 

TPAX(U.A.GEN.T). A valid combination implies that the following 

must be satisfied: 

(1) it's a delivery route (2) aircraft A can fly route R 

(3) start time >= unit's ALD (4) delivery date <= RDD + MAXLATE 

(5) positive movement requirement (6) Maximum aircraft payload for 

the route must be greater than 25% its full load capability; this is 

for cost efficiency. (7) Suppply of aircraft type >= 0 

DSETX(U,A,GEN,T)$DELR(GEN) = YES $( 1$DDR(GEN) AND 1$DRA(GEN,A) AND 

(ORD(T) GE ALD(U)) AND 

(ORD(T) LE (RDD(U) - ROUND(SUM(K$DRK(GEN,K}.CTIME(A,GEN,K))/24) 
+ MAXLATE)) AND 

(1$SUM((I.K)$(DRI(GEN,1)* DRK(GEN.K)), MOVEUE(U.LK)) OR 
ISSUM((I.K)$(DRI(GEN,1)*DRK(GEN,K)), MOVEPAX(U,1K)) ) 

AND (1$(MAXLOAD(A.GEN) GT (0,25 * ACLOAD(A."1")) )} 

AND SUM (TP$(ORD(TP) LE ORD(T)), SUPPLY(A,TP)) ) ; 


ee HK He He RF KH HK HF 


* DSETDY(A.GEN,T) is used to control the allowable combination 

* of aircraft, route. and start time for the decision variable 

* Y(A,GEN,T). An aJlowable combination is one in which the 

* aircraft-route is compatible and that the supply of aircraft 

* type to date is greater than 0). 

DSETY(A,GEN,T)$DELR(GEN) = YES $( ISDRR(GEN) AND 1$DRA(GEN,A) AND 
SUM (TP$(ORD(TP) LE ORD(T)), SUPPLY(A,TP)) ), 

*dcfine priccout parameters 


PARAMETER TOTX(U,A,GEN,T) price out value for x; 

TOTX(U,A,GEN,T)$(DSETX(U,A,GEN,T) AND DELR(GEN))=0.0; 

*computations for TOTX 

TOTX(U,A,GEN,T)$(DSETX(U,A,GEN,T) AND DELR(GEN))=TOTX(U,A,GEN,T)+ 

SUM(I$(DSETX(U,A,GEN,T) *DRI(GEN,I)), ACBALLM(A,L,T)); 

LOOP(TP, 

TOTX(U,A,GEN,T)$(DSETX(U,A,GEN,T) AND DELR(GEN))=TOT X(U,A,GEN,T)- 
SUM(K$(1$DSETX(U.A,GEN.T) AND I1$DRK(GEN.K) AND 
((ORD(T) + ROUND(CTIME(A,.GEN,K)/24)) EQ ORDETP))), 
ACBALK.M(A,K.TP)) ); 


TOTX(U.A.GEN.T)$(DSETX(U,A,GEN.T) AND DELR(GEN))=TOTX(U,A,GEN,T)- 
(FLTTIME(A,GEN)*ACURATE.M(A)); 
LOOP(TP, 
TOTX(U.A.GEN.T)$(DSETX(U,A,GEN,T) AND DELR(GEN) AND (ORD(T) LE ORD(TP)))= 
TOTX(U,A,GEN,T)+(K CONSUME(A,GEN,T,TP)*ACCONSUME. M(A,TP)) ); 





LOOP(TP. 
TOTX(U.A,.GEN.T)$(DSETX(U,A,GEN,T) AND DELR(GEN))=TOTX(U,A,GEN,T)+ 
SUM(AF$(DRI(GEN,AF)*DDR(GEN)), 
MOGREQ(AF,A)*(GTIME(A,'ONLD')/24)*MOGUTILITY.M(AF,T))+ 
SUM(AF$(I1$DRK(GEN,AF) AND 1$DDR(GEN) AND 
(ORD(T) + ROUND(CTIME(A,GEN,AF)/24) EQ ORD(TP))), 
MOGREQ(AF,A)*(GTIME(A,'OFFLD')/24)*MOGUTILITY.M(AF,TP))+ 
SUM(AF$(1$DRE(GEN,AF) AND 1$DDR(GEN) AND 
(ORD(T) + ROUND(CTIME(A,GEN,AF)/24) EQ ORD(TP))), 
MOGREQ(AF,A)*(GTIME(A,'ENR')/24)* MOGUTILITY.M(AF,TP)) ); 
*computation for TONSUE 


PARAMETER TOTTON(U,A,GEN,T); 

TOTTON(U.A.GEN,T)$(DSETX(U,A,GEN,T) AND DELR(GEN))=0.0; 
TOTTON(U.A,GEN.T)$(DSETX(U,A,GEN,T) ANF DELR(GEN))=TOTTON(U,A,GEN,T)+ 
SUM((I.K)$(ACSQFT(A) AND I$DSETX(U - .GEN,T) AND I$DRI(GEN,I) AND 

ISDRK(GEN,K) AND I$MOVEUE(U,I,K)), REQMTUE.M(U,I,K)); 


TOTTON(U,A,GEN,T)$(DSETX(U,A,GEN,T) AND DELR(GEN))=TOTTON(U,A,GEN,T)+ 
SUM((1,K)$(1SMOVEUE(U.I,K) AND ACCARGO(A,'OUT') AND 
ISDSETX(U.A,GEN,T) AND ISDRI(GEN,I) AND 1$DRK(GEN,K)), 
OUTREQMT.M(U,I.K)); 


TOTTON(U,A,GEN,T)$(DSETX(U,A,GEN,T) AND DELR(GEN))=TOTTON(U,A,GEN,T)+ 
SUM((I,K)$(1$MOVEUE(U,I,K) AND ACCARGO(A,'OVER') AND 
I1$DSETX(U,A,GEN.T) AND I$DRI(GEN,I) AND 1$DRK(GEN,K)), 
OVERREQMT.M(U,I,K)); 


*includes cost coefficient 
LOOP(K. 
TOTTON(U,A,GEN,T)$(ACSQFT(A) AND 1$DSETX(U,A,GEN,T) AND I$DELR(GEN) AND 
I$DRK(GEN,K))=TOTTON(U,A,GEN,T)- 
LATEPEN(U,'UE')*MAX(0,(ORD(T)+ROUND(CTIME(A,GEN,K)/24)-RDD(U))) ); 


*computations for variable 

PARAMETER TOTY(A,GEN,T); 

TOTY(A,GEN,T)$(DSETY(A,GEN,T) AND DELR(GEN))=0.0; 

LOOP(TP. 

TOTY(A,GEN,T)$(DSET Y(A,GEN,T) AND DELR(GEN))=TOTY(A,GEN,T)- 
SUM(I$(1$DRI(GEN,1) AND 
((ORD(T)+ROUND(CTIME(A,GEN,1)/24)) EQ ORD(TP))), ACBALI.M(A,],TP)) ); 


TOT Y(A,GEN,T)$(DSETY(A,GEN,T) AND DELR(GEN))=TOTY(A,GEN,T)+ 
SUM(K$(DSETY(A,GEN,T)*DRK(GEN,K)), ACBALK.M(A,K,T)); 


TOTY(A,GEN,T)$(DSETY(A,GEN,T) AND DELR(GEN))=TOTY(A,GEN,T)+ 
(FLTYIME(A,GEN)*ACURATE.M(A)); 
LOOP(TP, 
TOTY(A,GEN,T)$(DSETY(A,GEN,T) AND DELR(GEN) AND (ORD(T) LE 
ORD(TP)))=TOTY(A,GEN,T)+ 
(KCONSUME(A.GEN,T,TP)*ACCONSUME.M(A,TP)) ); 
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LOOP(TP. 

TOTY(A.GEN.T)$(DSETY(A,GEN,T) AND DELR(GEN))=TOTY(A,GEN,T)+ 
SUM(AF$(I$DSETY(A,GEN,T) AND 1$DRE(GEN,AF) AND I$DRR(GEN) AND 
($(ORD(T) + ROUND(CTIME(A,GEN,AF)/24) EQ ORD(TP))), 
MOGREQ(AF.A)*(GTIME(A.'ENR')/24)* MOGUTILITY.M(AF,TP)) ): 


*computations for variable TPAX 

PARAMETER TOTTPAX(U,A.GEN,T): 

TOTTPAX(U,A,GEN,T)$(DSETX(U.A,GEN,T) AND DELR(GEN))=0.0; 

TOTTPAX(U,A,GEN,T)$(DSETX(U,A,GEN,T) AND DELR(GEN))=TOTTPAX(U,A,GEN,T)+ 
SUM((I.K)$(1$DSETX(U,A,GEN,T) AND 1$DRI(GEN,I) AND 1$DRK(GEN,K) AND 
I$SMOVEPAX(U,I,K)), REQMTPAX,.M(U,I,K)); 

*include cost 

LOOP(K, 
TOTTPAX(U,A,GEN,T)$(DSETX(U,A,GEN,T)*DELR(GEN) * DRK(GEN,K))= 

TOTTPAX(U,A,GEN,T)- 
(LATEPEN(U.'PAX')*MAX(0,ORD(T)+ROUND(CTIME(A,GEN,K)/24)-RDD(U))) ): 


*warmup Starting 

X.M(U,A,GEN.T)$(DSETX(U.A,GEN.T) AND DELR(GEN))=1.0; 
Y.M(A.GEN,T)$(DSETY(A,GEN.T) AND DELR(GEN))=1.0; 
TONSUE.M(U,A,GEN,T)$(DSETX(U,A,GEN,T) AND DELR(GEN))=1.0; 
TPAX.M(U,A,GEN.T)$(DSETX(U,A,GEN,T) AND DELR(GEN))=1.0; 
ACWEIGHT.M(U,A,GEN,T)$(DSETX(U,A,GEN,T) AND DELR(GEN))=1.0; 
ACSPACE.M(U.A,GEN,T)$(DSETX(U,A,GEN,T) AND DELR(GEN))=1.0; 
ACPAX.M(U,A.GEN,T)$(DSETX(U,A,GEN,T) AND DELR(GEN))=1.0; 


*dcfine new route sct base to price out valucs 

LOOP((U.A.T), 

R(GEN)$(DSETX(U,A,GEN.T) * DELR(GEN))=YES$((TOTX(U,A,GEN,T) LT 0.0) OR 
(TOTTPAX(U.A.GEN,T) LT 0.0) OR (TOTTON(U,A,GEN,T) LT 0.0) OR (TOTY(A,GEN,T) LT 0.0)) ); 
DR(GEN)$DDR(GEN)=YES$R(GEN), 


RR(GEN)$DRR(GEN)=YES$R(GEN); 

RI(GEN,AF)$DRI(GEN, AF)=YES$R(GEN), 

RK(GEN, AF)$DRK(GEN, AF)=YES$R(GEN), 

RE(GEN, AF)$DRE(GEN, AF)=YES$R(GEN); 

ADJACENT(GEN, AF,AFP)$DADJACENT(GEN,AF,AFP)=YES$R(GEN), 
RA(GEN.A)$DRA(GEN, A)=YES$R(GEN). 

RESULT('NEWR’')=CARD(R). 

DISPLAY RESULT: 


SOLVE AIRLIFT USING RMIP MINIMIZING Z; 
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APPENDIX D. FLOW DIAGRAM OF ROUTE SELECTION MODEL 


START 


AGGREG. PAS (3) 






TPFDD. TXT (1) i 
ORIGIN.EXT (2) 
DESTINA.EXT (2) : | 










































TPEFDD. TXT 
ORIGIN.TXT 
DESTINA.TXT 


UNIT!. TXT 


DESTINA1 TXT 
UESQFT! TXT 
CARGOP] TXT 
ALDI. 1XT 
RDDI TXT 
ROUTELOG TXT (9 













HOFCENTROIDS 
K1 


AIRFIELD. SET 
DIST1. TXT 
#OFCENTROIDS (K) 


AGGREG. GMS (4) 





DESTSET. GMS 







ORIGSET. GMS 





DIESTINA 1, 1X1 
# OFCENTROIDS (KK) 










UNITI, DAT 
MOVING]. DAT 
MOVEPA11. DAT 
MOGCAP. TXT 
ENROUTE. TXT 








AIRFIELD. SET 
DIST1. TXT 

#OFCENTROIDS (K2) 
ORIGIN]. TXT 













MOGCAP. TXT 
ENROUTE. TXT 
MOVEMENT. TXT 
PAIRS. TXT 
DISTANCE. TXT 
MOGCAP. DAT 
ORIGIN. TXT 
DESTINA. TXT 
DISTANCE. DAT 
ROUTELOG. TXT 


SECAGGRE. PAS (5) 


ALDI. TXT 
RDDI. TXT 
MOVEMENT TXT 
CARGOPL.TXT 
UFESQFT1 TXT 
ORIGIN TXT 
DESTINA. TXT 
CONT... 





ENROUTE. TXT 
CREW. TXT 
AIRFIELD. TXT 
AIRCRAFT. TXT 
LAMBDA. TXT 





SECAGGRE. PAS 


UNITNAME. DAT 
MOVEUE. DAT 
MOVEPAX. DAT 
UESQFT. DAT 
CARGOP. DAT 
ALDT. DAT 
RDPT. DAT 
ROUTELOG. TXT 


ROUTESEL. PAS (6) 


PAIRS. TXT 
ENROUTE. TXT 
CREW. TXT 
AIRFIELD. TXT 
AIRCRAFT. TXT 
DISTANCE.TXT 
LAMBDA. TXT 


GAMSREG. DAT 
















ROUTELOG.TXT 


THROUHPUT. GMS (7) 









GAMSREG. DAT 
MOVEUE. DAT 
MOVEPAX. DAT 
DISTANCE. DAT 
MOGCAP. DAT 
UNITNAME. DAT 
UESQFT. DAT 
CARGOP. DAT 
ALD. DAT 
RDD. DAT 


COLUMN.GMS (8) 


GAMSCOL.DAT 
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Text file separated with tab and generated in EXCEL. 


Preselected origin-destination airfields which will be considered during the 
avgregation Ifthe user wants to include all the airfields that mentioned in TPFDD, 


then enter “No” for the including file question. 


Unit aggregation based on preselected origin destination airfields , minimum 


movement requirements and the difference between the RDD values 
Centroid model for origin and destination airfields written in GAMS. 
Unit aggregation based on centroid airfields. 

Route selection model. 

Optimization model. 


Sensitivity Analysis. 


Text file contains all the input and aggregation information based on this run. 
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